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ABSTRACT 

We have used previously published observations of the CO emission from the 
Antennae (NGC 4038/39) to study the detailed properties of the super giant 
molecular complexes with the goal of understanding the formation of young mas- 
sive star clusters. Over a mass range from 5 x 10^ to 9 x 10^ M©, the molecular 
complexes follow a power-law mass function with a slope of —1.4 ± 0.1, which is 
very similar to the slope seen at lower masses in molecular clouds and cloud cores 
in the Galaxy. Compared to the spiral galaxy M51, which has a similar surface 
density and total mass of molecular gas, the Antennae contain clouds that are an 
order of magnitude more massive. Many of the youngest star clusters lie in the 
gas-rich overlap region, where extinctions as high as ~ 100 mag imply that 
the clusters must lie in front of the gas. Young clusters found in other regions of 
the galaxies can be as far as 2 kpc from the nearest massive cloud, which suggests 
that either young clusters can form occasionally in clouds less massive than 5 x 10® 
Mq or these young clusters have already destroyed their parent molecular clouds. 
Combining data on the young clusters, thermal and nonthermal radio sources, 
and the molecular gas suggests that young massive clusters could have formed 
at a constant rate in the Antennae over the last 160 Myr and that sufficient 
gas exists to sustain this cluster formation rate well into the future. However, 
this conclusion requires that a very high fraction of the massive clusters that 
form initially in the Antennae do not survive as long as 100 Myr. Furthermore, 

^Department of Physics and Astronomy, McMaster University, Hamilton, Ontario L8S 4M1 Canada 
^Division of Physics, Mathematics, and Astronomy, Cahech 105-24, Pasadena CA 91125 U.S.A. 
^CEA/DSM/DAPNIA/Service d'Astrophysique, CE-Saclay, 91191 Gif sur Yvette Cedex, France 
^Astronomy Department, Cornell University, Ithaca NY 14853 U.S.A. 

''Chercheur Associe, Observatoire de Paris, LERMA, 61 Av. de I'Observatoire, F-75014 Paris, France 



- 2 - 



if most young massive clusters do survive for long periods, the Antennae must 
be experiencing a relatively short burst of cluster formation to prevent the final 
merger remnant from exceeding the observed specific frequency of star clusters 
in elliptical galaxies by a wide margin. Finally, we compare our data with two 
models for massive star cluster formation and conclude that the model where 
young massive star clusters form from dense cores within the observed super gi- 
ant molecular complexes is most consistent with our current understanding of 
this merging system. 

Subject headings: galaxies: individual (NGC 4038, NGC 4039, Arp 244) — galax- 
ies: ISM — ISM: molecules — stars: formation — radio lines: galaxies 

1. Introduction 

The discovery of very luminous young star clusters in merger remnant galaxies (Holtz- 
man et al. 1992; Whitmore et al. 1993; Whitmore & Schweizer 1995) has led to a dramatic 
shift in our understanding of star cluster formation. Some of these young star clusters may 
be massive enough (Zhang & Fall 1999; Ho & Filippenko 1996; Mengel et al. 2002) that they 
can be viewed as young counterparts to the ubiquitous globular clusters, which are found in 
nearly every extragalactic system from giant cD elliptical galaxies to dwarf irregular galaxies 
(Harris 2001). The discovery of these young massive clusters suggests that massive star 
cluster formation was not a process confined exclusively to the early universe (e.g. Cen 
(2001); Fall & Rees (1985)), but rather one that continues on to the present day. Thus, 
while attempting to understand the formation and evolution of young massive star clusters 
is interesting in its own right, such studies may also shed light on the process of globular 
cluster formation, which clearly occurred in a wide variety of galactic environments in the 
early universe. 

Since star clusters in the Milky Way are observed to form in dense clouds of molecular 
hydrogen gas, it is likely that these young massive star clusters formed in a similar way. 
However, the relatively large masses (> 10^ M0) estimated for these clusters are comparable 
to the masses of typical giant molecular clouds in nearby galaxies (Sanders, Scoville, & 
Solomon 1985). Since individual molecular clouds typically form stars with efficiencies of 
only a few percent (Evans & Lada 1991), these large cluster masses pose a challenge to our 
current understanding of star formation. In particular, if the star clusters are gravitationally 
bound, the star formation efficiency in the material from which they formed must have 
been closer to 50%. Such high star formation efficiency is typically observed only in smaller 
bound cores within individual giant molecular clouds (Lada et al. 1991a). These facts have 
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led to two quite different explanations for how these young massive star clusters formed. To 
explain the formation of globular clusters, Harris & Pudritz (1994) postulated the existence 
of "super giant molecular clouds" in the early universe. With masses up to 10^ Mq, but 
internal structures similar to local giant molecular clouds, these massive clouds could contain 
correspondingly massive cores, which could form massive (bound) star clusters. If similarly 
massive clouds exist in galaxies in the local universe, this might also be a viable method 
for forming young massive star clusters. An alternative explanation for young massive star 
clusters proposed by Schweizer et al. (1996) involves the collapse of a pre-existing giant 
molecular cloud to form a single massive star cluster. Such collapse could be triggered by 
overpressure in the diffuse interstellar medium produced during the collision of two massive 
galaxies (Jog & Solomon 1992), but would require a collapse of 2-3 orders of magnitude in 
linear size. Clearly, one way to shed some light on the processes which lead to the formation 
of young massive star clusters is to study the properties of the molecular gas in galaxies in 
which such clusters are found. 

Since the original discoveries of rich populations of young massive star clusters in nearby 
merging and merger-remnant galaxies, young massive clusters have been identified in a very 
wide range of galactic environments, from dwarf starburst galaxies (i.e. M82, O'Connell et 
al. (1995)) to spiral galaxies (i.e. M51, Larsen (2000)) to merging galaxies in all stages of 
evolution (see Whitmore (2003) for a complete list). The Antennae system (NGC 4038/39, 
Arp 244) has the largest number of young massive star clusters identified in a single galaxy 
to date (Whitmore & Schweizer 1995; Whitmore et al. 1999). In addition, the clusters in the 
Antennae are spread over a large area of the galaxy (~ 100 kpc^), which implies that these 
clusters have formed over much of the galactic disks. The Antennae have a higher surface 
density of young massive star clusters than older merger systems, although smaller galaxies 
such as M82 (O'Connell et al. 1995) have a larger absolute surface density of clusters confined 
to a relatively small region. This extremely rich population of clusters, combined with the 
proximity of the Antennae system (19 Mpc for Ho = 75 km s^^ Mpc^^), makes it the logical 
first choice for detailed studies of the structure and properties of the molecular interstellar 
medium with the aim of understanding the formation of young massive star clusters. 

The first detection of CO in the Antennae by Sanders & Mirabel (1985) suggested that 
the system had a relatively high ratio of far-infrared to CO luminosity. More recent com- 
plete mapping of the two galactic disks by Gao et al. (2001) revealed a much larger reservoir 
of molecular gas, which led to the suggestion that this system could possibly evolve into 
an ultraluminous starburst as the merging becomes more advanced. Early interferometric 
CO observations of the Antennae system were presented by Stanford et al. (1990). These 
observations revealed that the most massive concentration of molecular gas lay outside the 
two nuclei in a region which Stanford et al. named "the overlap region". In a recent paper 
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( Wilson et al. 2000, hereafter Paper I), we presented interferometric CO observations of 
the Antennae with improved resolution and sensitivity and covering a larger fraction of the 
galactic disks. These data show the presence of massive (> 10*^ Mq), gravitationally bound 
gas clouds, which we term super giant molecular complexes (Paper 1). A detailed compari- 
son of the CO data with high-resolution mid-infrared data from the Infrared Astronomical 
Satellite (Vigroux et al. 1996; Mirabel et al. 1998) suggests that the extremely bright mid- 
infrared source seen in the overlap region could be produced by star formation induced by 
the collision of two of these massive complexes (however, see Liang et al. (2001) for evidence 
that the clouds may not be colliding) . An alternative explanation is that this region contains 
extremely young (< 1 Myr) sites of star formation (Paper I). 

In this paper, we use the data set presented in Paper I for a detailed study of the 
population of molecular complexes in the Antennae, such as their mass function, their spatial 
overlap with the young massive star clusters, and a comparison with theories of star cluster 
formation. In §2 we describe the observations and data reduction, as well as a brief discussion 
of how the clouds were identified from the data cube. In §3 we present the cloud mass 
spectrum and compare our results to observations of massive cloud complexes in other nearby 
galaxies. We also discuss what can be learned from the relative locations of the young massive 
star clusters and the super giant molecular complexes. In §4, we make a detailed comparison 
of the Antennae with the complexes and star clusters observed in M51 (Rand Sz Kulkarni 
1990; Larsen 2000; Scoville et al. 2001). In §5 we discuss the recent star formation history 
in the Antennae and in §6 we compare the properties of the molecular complexes with two 
models for massive star cluster formation. The paper is summarized in §7. 



2. Analysis 

2.1. Observations and Data Reduction 

The CO observations of the Antennae originally presented in Paper I were obtained with 
the Caltech Millimeter Array between 1998 March and 1999 February. Three overlapping 
field centers were chosen to cover the majority of the young massive star clusters identified 
by Whitmore & Schweizer (1995). The coordinates of the three field centers were specified 
in 1950 coordinates as (11:59:20.3, -18:35:31) (Field 1), (11:59:20.3, -18:36:11) (Field 2) and 
(11:59:17.5, -18:35:21) (Field 3), which gives an offset between adjacent field centers of 40". 
In J2000 coordinates, the equivalent field centers are (12:01:54.09, -18:52:13.25), (12:01:54.09, 
-18:52:53.25), and (12:01:51.28, -18:52:03.25). Each field was observed for the equivalent of 
one complete transit in each of four configurations: C (compact), L (low), E (equatorial) 
and H (high resolution). The baselines for these observations ranged from 15 m to 250 m. 
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The observations of Fields 1 and 2 were interleaved in each observing track, which results in 
very similar uv coverage for these two fields. Field 3 (which was added later) was observed 
in separate tracks. These CO data represent approximately 90 hours observing time with 
the array. 

The spectrometer was configured to give 2 MHz (5.2 km s~^) resolution and was tuned 
to a central velocity of 1547 km s""*-; the resulting velocity range was from 1256 km s~^ 
to 1838 km s^-*^. The system temperatures obtained during these observations range from 
1000 to 2000 K (single sideband), with an average value of 1500 K. (Note that these system 
temperatures arc substantially higher than those typically obtained by the array due to the 
low declination of the source.) The data calibration was performed using the mma package 
(Scoville et al. 1993). The bright quasar 3C273 was observed for flux cahbration; its flux was 
derived from observations of either Uranus or Neptune taken from other observing sessions in 
the array database. The flux of 3C273 changed slowly during the 12 months over which these 
data were obtained. The following values were derived for the flux of 3C273: 21 Jy (1998 
March-May); 22 Jy (1998 September - November 1); 19.5 Jy (1998 November 1 - December 
20); 17 Jy (1998 December 20 - 1999 February 16). The absolute calibration uncertainty 
is estimated to be 20%. The gain and passband calibrator for these observations was the 
bright quasar 3C279. Only observations for which the coherence on 3C279 was measured to 
be > 50% were included in the flnal analysis. 

The data were mapped using the Miriad analysis package (Sault, Teuben, & Wright 
1995). The uv data were first clipped to have a maximum amplitude of 14 Jy. The three 
fields were inverted together with robust weighting to make a single mosaic dirty map. The 
rms noise in this dirty map in line-free channels was 0.055 Jy beam~^ and the synthesized 
beam was 3.15 x 4.91" (or 310 x 480 pc at a distance of 19 Mpc). The map was cleaned using 
the task MOSSDI with a clean cutoff of 0.11 Jy beam~^ (2c). We used a single clean box 
made by summing three 60 x 60" clean boxes centered on each of the three field centers. As 
a final step, individual maps and the data cube were corrected for the fall-off in sensitivity 
due to the primary beam of the 10.4 m antennas. We used the task MOSSEN to create a 
map of the gain due to the primary beam and then divided individual maps by that gain 
file, while also masking regions of the maps where the gain map was smaller than 0.5. 



2.2. Identifying Clouds from the Data Cube 

Individual clouds were identified from the data cube using the automatic clump iden- 
tification algorithm CLFIND (Williams, de Geus, & Blitz 1994). This algorithm searches 
for peaks of emission within a contour map of the data, which it then follows down to lower 
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intensity levels, and has the advantage of not assuming any specific clump profile (such as 
a gaussian profile). The main free parameter in the algorithm is the contour level, which 
defines our cloud detection threshold. The contour level is specified relative to the noise in 
the map, and is usually set to 2a. Only emission brighter than one contour level can be 
included in a cloud, and a cloud must have at least one pixel that is twice the minimum 
contour level in order to be found. Thus, the initial cloud identification was run on the data 
cube before correction for primary beam attenuation. We experimented with three different 
contour levels to see how the contour level affected the properties of the clouds; we used 
contour levels of 0.10, 0.11 (2(7), and 0.12 Jy beam~^. Unhke Paper I, individual clouds were 
not inspected and merged by hand where they appeared to overlap; thus, the properties of 
the largest clouds will be somewhat different from the values in Paper I. 

The properties of the individual clouds such as position, flux, and velocity were mea- 
sured using the program CLSTATS (WiUiams et al 1994). This program combines the clump 
assignment cube produced by CLFIND with the original data cube to calculate the parame- 
ters of the individual clouds from the pixels that have been assigned to that cloud. We used 
the gain-corrected data cube in running CLSTATS to produce fluxes that were corrected for 
the primary beam. The properties of the clouds measured with a contour level of 0.11 Jy 
beam^^ arc given in Appendix A. 

To estimate the sensitivity of the flux of an individual cloud (and hence its mass) to the 
chosen contour level, wc took a closer look at the largest clouds in our sample, those with 
masses estimated to be larger than 10^ Mq. Nine large clouds idcntiflcd with a contour level 
of 0.11 Jy beam^^ were cross-identified with the lists from the other two contour levels, and 
the masses measured for each of the three identification runs were compared. The average 
dispersion in the mass measurements determined from the three different identification runs 
was 25% while the dispersion for an individual clouds ranged from a low of to a high of 
55%. Thus, we estimate the random uncertainty on the mass estimate for any individual 
cloud in our sample to be 25%. 



2.3. Moment Maps and Total Flux 

The integrated intensity map presented in Paper I was made using the command 
CLPLOT to sum together the emission in the data cube that corresponded to molecular 
complexes identified by CLFIND. Only molecular complexes with emission in at least three 
velocity channels were included in the integrated intensity map. Thus, this image (shown 
again in Figure la) differs from a standard zeroth moment map because very small clouds and 
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isolated, weak emission regions have not been included^. An integrated intensity map (not 
shown) produced in the same way but including all the complexes identified by CLFIND, no 
matter how weak, is very similar to the positive contours in the zeroth moment map made 
with a 2(7 intensity cutofi^ (Figure Ic). Figure 2 shows the first moment map of the molecular 
complexes in the Antennae. The original data cube was filtered by keeping only emission 
that CLFIND identified with a molecular complex; the first moment map was produced from 
this filtered data cube with no additional flux cutoff. Thus, this flrst moment map includes 
emission down to a level of 2a, but is somewhat cleaner than a flrst moment map made with 
a 2(7 cutoff from the original data cube because isolated emission regions extending over 
fewer than four pixels have not been included. 

The zeroth moment map shows clear negative bowls (which are not shown in Figure 1 for 
clarity) around the bright emission regions, which are a sign that we are not detecting all the 
flux in the region with the interferometer. These negative bowls make it difficult to measure 
the total flux, so the total fiux was measured instead from the integrated intensity image from 
CLPLOT which included the very small clouds. The integrated intensity measured in this 
way is 930 Jy km s~^, compared to 910 Jy km s^^ measured in Paper I. This flux is probably 
a slight underestimate of the total flux in the map, since CLFIND sometimes misses weak 
emission at the edges of clouds. To correct for this effect, we calculated the difference between 
the CLFIND map and the zeroth moment map, and measured an additional flux of 125 Jy 
km s~^ in the two nuclei and the overlap region where signiflcant positive signal was seen. 
Thus, the flnal estimate of the total integrated intensity detected with the interferometer 
is 1055 Jy km s~^. Adopting a CO-to-H2 conversion factor of 3 x 10^° H2 cm~^ (K km 
s~^)^^ (Strong et al. 1988) and including a factor of 1.36 to account for heavy elements, this 
intensity corresponds to 6.1 x 10^ M© of molecular gas. In Appendix B, we estimate the 
CO-to-Il2 conversion factor for the Antennae from a few large clouds and flnd that it agrees 
with the Galactic value within the uncertainties. 

Gao et al. (2001) have made a complete single-dish map of the Antennae system. They 
measure a total flux of 3172 Jy km s^^, which, for the same distance and conversion factor 
adopted above, corresponds to a total mass of 1.8 x 10^'^ Mq. This flux is signiflcantly larger 
than the flux detected in our interferometric map. However, their map, which contains 73 
spectra at 28" spacing, covers a substantially larger area than does our map. To estimate 
the single dish flux within our smaller map, we summed the spectral lines from fourteen 



^The version of CLPLOT that was used in this analysis had a bug in that the last cloud of a list of 20 
was not included in the plots. This bug resulted in one small cloud not being included near the southern 
end of the west ring and some small changes to the contours around the two nuclei. The corrected version 
of the plot is shown in Figure lb. 
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spectra which fell within the clean box used in processing the interferometer map. The total 
flux from these fourteen spectra is 1654 Jy km for a mass of 9.6 x 10^ Mq. Thus, the 
interferometer map appears to have detected 65% of the total flux seen over the same area 
in the fully-sampled single-dish map. 

3. Molecular Complexes in the Antennae and Other Galaxies 

3.1. The Mass Spectrum of the Moleculeir Complexes 

Wc measured the mass spectrum for the molecular complexes in the Antennae using 
the three different contouring levels discussed in §2 to see whether the slope and upper-mass 
cutoff to the mass spectrum were heavily dependent on the contour used to identify the 
clouds. The total number of complexes identified was 114 for a contour level of 0.10 Jy 
beam~^, 86 for a contour level of 0.11 Jy beam~^, and 49 for a contour level of 0.12 Jy 
beam~^. For each set of complexes, we calculated the number of the clouds in logarithmic 
mass bins with a spacing of 0.2; two different centers for the mass bins were also used, with 
the first having a bin centered at log(M) = 7.0 and the second having a bin centered at 
log(M) = 7.1. The number of clouds in each mass bin was then divided by the central 
mass of that bin to obtain a true differential mass function, dN / dM . The differential mass 
functions for the three different identification runs are shown in Figure 3. 

We calculated the detection and completeness limits for our data using a modified 
version of the formulae given in Hcycr ct al. (2001). For these inter fcromctric data processed 
using CLFIND, a cloud must contain at least four pixels, one of which exceeds twice the 
specified contour level, in order to be detected. The detection limit in Jy km s^^ is then 
given by 

S'S^ = 5A,5Ay/17.55 

where is the contour level in Jy beam^^, AV^ is the velocity width of a single channel (5.2 
km s~^) and the factor of 17.55 is the area of the synthesized beam in pixels. For A^" = 0.11 
Jy beam~^, the detection limit is Sqq' — 0-16 Jy km s~^. The 5(7 completeness hmit is given 
by 

^co = ^co + 5o-('S'co) 

where 

(y{Sco) = aAV^NpNjl7.55 

and a is the rms noise in the map Np is the minimum number of pixels required for a 
detection and is the minimum number of velocity channels required for a detection. For 
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(7 = 0.055 Jy beam ^, Np — 4, and Nc — 1, we get a{Sco) — 0.14 Jy km s ^, which gives a 
5(7 completeness hmit of 0.86 Jy km s~^ or a mass hmit of 5.0 x 10^ Mq. 

We made least squares fits to the six mass functions including \fN uncertainties and 
ignoring bins with log(M) < 6.7. The slope and uncertainty in the slope are included in 
Figure 3. The average slope across the three different identification runs is a = —1.4 ± 0.1. 
This slope refers to a mass range from 5 X 10^ Mq to ~ 10*^ Mq and is significantly shallower 
than the canonical stellar initial mass function derived by Salpeter (1955) (-2.35 on the 
same scale). However, this slope of -1.4 is very similar to the slope of -1.5 obtained for the 
mass function of Giant Molecular Clouds in the Milky Way (Sanders, Scoville, & Solomon 
1985; Solomon et al. 1987). It is somewhat shallower than the slope of -1.8 found over a 
mass range of 1000 M© to 10^ for molecular regions in the outer Galaxy (Heyer et al. 
2001). However, Heyer et al. (2001) point out that many of the low-luminosity objects in 
their survey are not self-gravitating, in which case the mass function may be shallower than 
the luminosity function. In addition, the relatively poor spatial resolution of our Antennae 
data means that blending of individual clouds into a single feature could be an important 
effect, which would tend to make our observed mass function shallower than the true mass 
function. Kramer et al. (1998) also obtained similar slopes ranging from -1.6 to -1.8 from 
GO observations of clumps inside seven molecular clouds. The clump masses in their sample 
ranged from as low as 10~^ Mq to as high as 10^ Mq and overlap at the high-mass end with 
the mass range covered by Heyer et al. (2001). [However, steeper clump mass functions, 
consistent with the stellar initial mass function over a fairly narrow mass range from ~ 0.5 
Mq to ~ 10 Mq, have been derived recently for several nearby clouds using submiUimeter 
continuum data (Testi & Sargent 1998; Mottc ct al. 1998; Johnstone et al. 2000; Motte et 
al. 2001).] Thus, by combining our new study of the Antennae with previous work on the 
Milky Way, we see that the mass function of structures in the molecular interstellar medium 
in galaxies has a fairly constant slope of ~ —1.5 ± 0.2 that extends over at least 8-9 orders 
of magnitude in mass, from 1-10 Mq up to 10^ Mq. It is worth noting that these power-law 
mass spectra are consistent with a model where clouds grow by agglomeration from smaller 
objects (Kwan 1979). 



3.2. Previous Observations of Moleculeir Complexes in Other Galetxies 

We have searched the literature for GO observations of galaxies with spatial resolution of 
better than 1 kpc to see what is known about massive molecular complexes in other galaxies 
(Table 1). Objects similar to the most massive complexes in the Antennae have been found 
in only three other galaxies: NGG 1068 (Planesas et al. 1991), Arp 220 (Sakamoto 1996), 
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and Arp 299 (Casoli et al. 1999). Adopting a distance of 14.4 Mpc for NGC 1068 (TuUy 
& Fisher 1988), the most massive complex is 4.5 x 10^ (Planesas et al. 1991) and the 
average molecular gas surface density in the inner arcminute is roughly 100 Mq pc~^, both 
comparable to what is seen in the Antennae. In Arp 220 and Arp 299, the massive gas 
concentrations are located in the galactic nuclei, which probably makes them more similar 
to the gas located in the nucleus of NGC 4038 than the molecular complexes in the disks 
and overlap region. In the remaining normal spiral galaxies for which high resolution CO 
observations have been made, the most massive complex ranges from as small as 1.6 x 10^ 
Mq in M83 (Rand et al. 1999) to as large as 1.3 x 10^ Mq in MlOO (Rand 1995). For the five 
galaxies in Table 1 with distances between 9 and 19 Mpc, the angular resolution of the data 
sets varies by only a factor of two, while the mass of the most massive complex varies by a 
factor of almost twenty. This result suggests that the large complexes seen in the Antennae 
and in NGC 1068 are not the result of cloud blending due to the relatively large distances 
of these two galaxies. 

In only two galaxies, M51 (Rand & Kulkarni 1990) and NGC 1068 (Planesas et al. 
1991), are sufficiently large numbers of objects detected that it is possible to look at the 
mass function of the molecular complexes. However, the limited mass range of the molecular 
complexes in M51 (1 — 6 x 10^ M©) makes it impossible to determine the slope with any 
degree of accuracy. The 38 complexes detected in NGC 1068 range from 2 x 10^ Mq to 7 x 10^ 
Mq. Fitting the entire mass range in the same manner as was done for the Antennae data 
gives a slope of —1.3 ± 0.2 for the differential mass function, which is quite similar to that 
obtained in the Antennae. From this comparison with previous studies of other galaxies, it 
is clear that the population of molecular complexes identified in the Antennae is both the 
largest sample and covers the largest dynamic range in mass. It would be interesting to carry 
out deeper observations of Arp 220 and Arp 299 to see if the mass functions of the molecular 
complexes in these more advanced interacting systems are similar to those measured for the 
Antennae and NGC 1068. 



3.3. The Spatial Relationship Between Young Massive Stair Clusters and 
Super Giant Moleculeir Complexes in the Antennae 

Whitmore et al. (1999) have published a large sample of young massive star clusters 
observed with WFPC2 on the Hubble Space Telescope. They find that the overlap region 
contains primarily clusters with ages < 5 Myr. This region contains roughly half the total 
amount of molecular gas detected in our map and has correspondingly very high extinctions, 
with an average Ay of 96 mag towards the three brightest CO peaks in this region. Thus, any 
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clusters identified in this part of the overlap region must lie on the near side of the molecular 
complexes. The western loop of gas and young stars to the west of NGC 4038 contains 
clusters with ages primarily in the range of 5-10 Myr. In this region, the CO emission occurs 
primarily between regions of many young clusters (Figure 4), and the typical extinction 
towards a cloud in this region is ~ 13. A very interesting region from the point of view 
of cluster formation is the north-eastern portion of NGC 4038. Whitmore et al. (1999) find 
that this region contains many young massive star clusters, with two-thirds of the clusters 
having ages less than 30 Myr and one-third having ages around 100 Myr. The presence of the 
younger star cluster population is interesting given the total lack of detected CO emission in 
this part of our map. Our 5a detection limit of 5 x 10^ Mq (§3) means that only relatively 
normal giant molecular clouds can be present in this region of the galaxy. 

Zhang, Fall, & Whitmore (2001) have compared the locations of three age groups of star 
clusters with data at a variety of other wavelengths. They find that clusters with ages < 5 
Myr (the reddest clusters) are more associated with wavelengths longer than mid-infrared, 
while clusters with ages > 10 Myr are more closely associated with far-ultraviolet and X-ray 
emission. The 2-point correlation function for the red clusters is a power-law up to a radius 
of 0.74 kpc. The maximum size of the power-law portion of the correlation function is very 
similar to the sizes of the super giant molecular complexes, which suggests that an upper 
hmit to the correlation function could correspond to a set of clusters formed inside a single 
such complex. Zhang et al. (2001) estimate a total star formation rate from Ha emission 
(corrected for 1.8 magnitudes of extinction and assuming a Salpeter initial mass function 
from 0.1 to 100 Mq) of 20 Mq yr~^. This star formation rate is a factor of two larger than 
the star formation rate of 10 Mq yr"-*^ which can be derived from the far-infrared luminosity 
following Kennicutt (1998). 

We compare the locations of the young, red star clusters discussed in Zhang et al. (2001) 
with the distribution of the molecular gas in Figure 4. Roughly 40% of these young, obscured 
clusters lie within the region of strong CO emission in the overlap region, which occupies 
only 10% of the total area mapped in CO. Thus, we confirm the conclusions of Zhang et 
al. (2001) that the youngest clusters are significantly correlated with regions containing 
molecular gas. However, some of these youngest star clusters lie a significant distance from 
strong CO emission, up to 1-2 kpc in regions to the north and east of the galaxies. This 
result is also true for the older Bl and B2 cluster groups [compare Figure 4 with Figure 2 of 
(Zhang et al. 2001)]. Since our CO map does not contain all the CO flux of these galaxies, it 
is possible that signiflcant molecular gas exists near these otherwise isolated young clusters. 
However, this molecular gas must be in rather low-mass clouds to have avoided detection 
in our map; in particular, any remaining clouds associated with these clusters are likely to 
have masses less than 5 x 10^ Mq, our 5a detection limit for clouds. Comparing the CO and 
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cluster distributions thus suggests that either relatively low-mass molecular clouds can form 
the occasional massive star cluster, or that massive star clusters can destroy their parent 
clouds on a rather rapid time scale, or both. 

4. A Case Study Comparison with M51 

To illustrate better the unusual aspects of the molecular interstellar medium in the 
Antennae, we compare the total gas mass and the properties of the molecular complexes and 
young massive clusters to the grand-design spiral galaxy M51. (Unfortunately there seem 
to be no surveys for young massive star clusters in NGC 1068, which would otherwise be 
another interesting galaxy to compare to the Antennae.) Rand & Kulkarni (1990) mapped 
a 4 X 5' region of the inner disk of M51 with the thrcc-clcmcnt Caltech Millimeter Array and 
detected 26 giant molecular associations with masses ranging from 1 x 10^ Mq to 5 x 10^ 
Mq. The beam size (9 x 7" or 370 x 290 pc at a distance of 8.6 Mpc, (Feldmeier, CiarduUo, 
& Jacoby 1997; Ferrarese et al. 2000)) is a very good match to our study of the Antennae. 
To match the total area surveyed as well, we restrict our analysis to clouds within a 5 kpc 
radius of the nucleus of M51. Correcting for the larger distance to the Antennae, our study 
is 1-2 times more sensitive to a complex of a given mass than is the M51 study of Rand & 
Kulkarni (1990). To correct for this difference in mass sensitivity, we further restrict our 
comparison to complexes that have masses larger than 2 x 10'^ Mq. More recent work by 
Aalto et al. (1999) provides better resolution and sensitivity for M51, but we prefer to use 
the older study of Rand & Kulkarni (1990) because it provides a better match in spatial 
resolution and mass sensitivity to our Antennae data. 

M51 contains 13 molecular complexes with masses greater than 2 x 10^ Mq inside an 
area of 79 kpc^. In contrast, the Antennae contain 43 complexes with masses greater than 
2 X 10^ Mq inside the same area. Thus, the Antennae have a three times higher surface 
density of massive molecular complexes than does M51. (However, correction for the poorly 
known inchnation of the disks of NGC 4038 and NGC 4039 might reduce this ratio by up to a 
factor of two.) In addition, the most massive complex seen in M51 has a mass of 5 x 10^ Mq, 
while the most massive complex in the Antennae (ignoring the gas in the nucleus of NGC 
4038) has a mass of 6 x 10* Mq, 12 times more massive than what is seen in M51. Thus, 
while the Antennae have a somewhat higher surface density of massive molecular complexes 
than docs M51, the most striking aspect of the molecular complexes in the Antennae is that 
their mass function extends to much larger masses than are seen in M51. If the molecular 
complexes in both galaxies follow a similar power-law mass function, then fitting the data 
with truncated power-laws following the prescription of Scoville et al. (2001) shows that the 
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M51 mass function is significantly truncated. In particular, given that 13 complexes are 
observed with masses greater than 2 x 10^ Mq, we would expect to see 4 complexes with 
masses greater than 6 x 10^ Mq. Thus, it is probably not small number statistics that is 
imposing the upper mass limit in M51. 

It is important to compare the total amount of molecular gas in the Antennae and M51 
to see whether the presence of very massive complexes can be attributed to a larger total 
reservoir of molecular gas to form them. Within the region of our study, Gao et al. (2001) 
measure a total mass of 9.6 x 10^ Mq (§2). This mass corresponds to a surface density of 
120 Mq pc~^ (again with no correction for inclination, which could reduce this value by up 
to a factor of two). Garcia-Burillo et al. (1993) made a large map of M51 using the IRAM 
30 m telescope; within the central 79 kpc^ and using a conversion factor that is three times 
smaller than the value adopted here for the Antennae, their measurements correspond to a 
total mass of molecular gas of 6 x 10^ M© or a surface density of 90 Mq pc~^. Thus, within 
uncertainties due to inclination in the Antennae and the CO-to-H2 conversion in M51, the 
total molecular gas surface density appears comparable in the two galaxies and certainly not 
different enough to account for the formation of extremely massive molecular complexes in 
the Antennae. 

One important difference between the two galaxies is their velocity fields. Figure 2 shows 
that the Antennae have a reasonably regular velocity field in the north-western arc, which 
may represent a fairly unperturbed piece of the disk of NGC 4038. However, the velocity 
field in the overlap region appears extremely disturbed; instead of a smooth progression in 
velocity, there are two large regions with fairly constant velocities separated by a boundary 
with a steep velocity gradient. In contrast to the disturbed velocity field of the Antennae 
system, the velocity field in M51 shows evidence of streaming motions due to the presence of 
spiral density waves superimposed on a very regular rotation pattern (Aalto et al. 1999). This 
marked difference in the velocity fields in M51 and the Antennae (which is not unexpected 
given the strong interaction in the Antennae) suggests a possible explanation for the presence 
of very massive molecular complexes in the Antennae. In a regularly rotating disk, molecular 
clouds that contain enough mass to be gravitationally bound may be subject to disruption 
by the effects of tides and shear (Rand 1993). Thus, it is possible that the formation of 
very massive gas clouds is suppressed in M51 due to high shear in its differentially rotating 
disk. The effect of galaxy mergers like the Antennae on the shear in the galactic disks is 
not well understood, althought it seems possible that some regions undergo increased shear 
while other regions may experienced reduced shear. It is interesting that for NGC 1068 (the 
other nearby spiral galaxy with massive gas clouds), Schinncrer et al. (2000) have suggested 
that the gaseous spiral arms may lie at the Inner Lindblad Resonance of a much larger (17 
kpc) bar. Thus, the massive gas clouds in NGC 1068 may also have formed in a region of 
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r educed shear. 

Surveys for young massive star clusters in M51 have been carried out recently by Larsen 
(2000) and Lamers et al. (2002). Using ground-based imaging, with HST archival confir- 
mation of 10 clusters, Larsen (2000) identified a population of 69 candidate young clusters 
with ages less than 500 Myr. This sample contains 40 clusters with My < — 10 mag and 
4 clusters with My < -12 mag. Lamers et al. (2002) used HST WFPC2 imaging of the 
nuclear region of M51 to identify a sample of 30 point-like sources. However, they concluded 
that these sources were most likely to be single stars or very small (M < few 100 Mq) star 
clusters, rather than massive star clusters. At the extremely young end of the scale, Scoville 
et al. (2001) concluded from their Paa and Ha surveys that the most massive star cluster in 
M51 that is ionizing an HII region has a mass < 5000 Mq. Larsen (2000) notes that M51 is 
rich in young massive star clusters compared to an average spiral galaxy, but that it is not 
unique in its rich star cluster population. 

In comparison, in just a single 10 kpc^ region of the Antennae, Whitmore et al. (1999) 
find a total of seven clusters with My < —12 mag. Although this region of the Antennae 
is particularly rich in clusters, it seems likely that the Antennae as a whole contain an 
order of magnitude more of the most luminous clusters than does M51. A total of 14,000 
point sources have been identified in the Antennae (Whitmore et al. 1999), of which at least 
40% are contaminating stars and at least 800 are definitely star clusters. Zhang & Fall 
(1999) identify 2000 star cluster candidates from this sample which lie within 0.3 mag of 
cluster evolutionary tracks and have Af,u,o < —9 mag. Compared to the ~ 70 star clusters 
identified in M51 by Larsen (2000), these observations suggest that the Antennae contain 
at least 30 times as many star clusters as M51. This difference is even more striking when 
we consider the average star formation rates in the two galaxies. The star formation rate 
is roughly proportional to the far infrared luminosity (i.e. Kennicutt (1998)); M51 has a 
far-infrared luminosity of 1.8 x 10^° Lq (Rice et al. 1988), while the far-infrared luminosity 
in the Antennae is 5.6 x 10^'^ Lq (Gao et al. 2001), roughly 3 times larger than M51. Thus, 
in comparison to their star formation rates, the Antennae contain an order of magnitude 
more young massive star clusters than does M51. We suggest in §6 that the formation of 
young massive clusters in the Antennae is enhanced by its ability to form very massive gas 
clouds. 



- 15 - 



5. Is the Antennae Undergoing a Short-Lived Steirburst? 

5.1. Recent Star Cluster Formation Rates from Optical and Radio Data 

We can derive clues to the recent star cluster formation rate in the Antennae by ex- 
amining the number of massive star clusters as a function of age. We can also use radio 
continuum observations to estimate the number of very young, embedded star clusters and 
older supernova remnants (Neff & Ulvestad 2000). The mass function analysis of Zhang & 
Fall (1999) shows that there are roughly 300 young star clusters more massive than 10^ Mq 
in the Antennae. Of these clusters, roughly 100 fall in the young age range of 2.5-6.3 Myr 
and a similar number fall in the age range of 25-160 Myr. There arc two simple explanations 
for the similar number of star clusters observed in these very different size age bins. If the 
observed number of star clusters reflects the average cluster formation rate over each time 
period, then the cluster formation rate in the last 6 Myr must have been 20 times higher 
than the cluster formation rate over the previous 150 Myr. However, this picture implies 
a very rapid change in the cluster formation rate compared to the current crossing time of 
the system (~ 100 Myr, Mihos, Bothun, & Richstone (1993)). Alternatively, if the cluster 
formation rate has been constant over the last 150 Myr, then only roughly one out of twenty 
of the younger clusters can survive in the long term (Zhang & Fall 1999; Whitmore 2002). 
Mengel et al. (2002) have measured dynamical masses for five bright clusters in the Anten- 
nae, which suggest that these clusters are gravitationally bound and hence likely to survive. 
However, we have no information about the likelihood of survival for the vast majority of the 
star clusters in the Antennae. A counter example of a luminous star cluster that seems to 
have insufficient low- mass stars to survive can be found in M82 (Smith & Gallagher 2001). 

A recent high-resolution radio survey shows that the Antennae contain 13 sources with 
obviously thermal spectral energy distributions which have likely ages in the range of 1-3 
Myr (Neff & Ulvestad 2000). These sources contain an equivalent mass in stars between 0.1 
and 60 Mq (assuming a Salpeter initial mass function) of 2 x 10^ to 1 x 10^ M© Roughly 
accounting for sources that have been missed because of large uncertainties on their spectral 
index as well as slightly less massive sources, we estimate that the number of optically thin 
thermal sources with masses greater than 10^ Mq is perhaps 30. In comparison, there are 
roughly 100 star clusters more massive than 10^ Mq with ages < 6 Myr (Zhang & Fall 1999). 
Thus, the number of optically visible star clusters and the number of optically thin thermal 
sources appear to be consistent with a constant star cluster formation rate in the Antennae 
over the last 6 Myr. 



''These masses can be obtained by scaling the mass in 05 stars by a factor of 45; assumes ionizing 
luminosityLj cx mP'^, Lj(05) = 4.7 x 10*® erg s~^, and M(05) = 60 M©; Ulvestad, private communication. 
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These numbers are also consistent with the result from Whitmore & Zhang (2002), who 
find that essentially all the bright thermal radio sources can be identified with optically visible 
star clusters. This good match between the radio and optical sources is rather puzzling in the 
overlap region, where the visual extinction can reach as high as 100 mag (§3.3). However, a 
closer examination of the sources in the overlap region reveals that all sources except WS-80 
lie well away from the brightest CO peaks. WS-80 itself is coincident with SGMC 7, which 
has an average visual extinction of 25 mag. Whitmore & Zhang (2002) estimate a visual 
extinction for WS-80 of 7.6 mag, which is consistent with it lying fairly close to the center 
of this somewhat lower mass molecular cloud. 

Turning now to the supernova rate, which traces star formation over slightly longer 
time periods, Neff & Ulvestad (2000) used the observed flux in non-thermal radio sources 
to estimate a current supernova rate of 0.1-0.2 yr~^. This rate is a factor of 10-20 larger 
than the rate of 0.01 yr~^ they estimate from the number of 05 stars and the lifetime of 
these stars. This result led them to postulate that the Antennae suffered a very sharp burst 
of star formation 3-4 Myr ago that lasted only a few hundred thousand years. However, 
the 60 Mq 05 stars are not the only stars that will produce supernovae; conservatively, all 
stars with masses greater than 20 Mq will produce a supernova. There will be many more 
of these lower mass O stars but they will have significantly longer main sequence lifetimes. 
However, if the star formation rate has been constant over the longer hfetimes of these 
stars, then the numbers of O stars will approach a steady state, with the number of stars 
at any mass proportional to the product of the initial mass function and the lifetime of the 
stars. This effect will increase the supernova rate over that calculated by Neff & Ulvestad 
(2000). For example, if the stellar lifetime were to vary as t ~ m~^'^^, then the number of 
supernovae produced by 20-50 Mq stars would be 35 times the number produced by 50-60 
Mq stars. Thus, depending on the exact relationship between mass and lifetime and the 
lower mass limit for producing a supernovae, it appears that the observed supernova rate in 
the Antennae is roughly consistent with a constant star formation rate over the last 10-20 
Myr. 

Given the results from the comparison with the radio data, it seems reasonable to assume 
that the star cluster formation rate over the last 160 Myr was roughly constant, but that 
only a small fraction of the clusters formed 160 Myr ago survive to the present time (Fall & 
Zhang 2001). This scenario would involve forming ~ 2000 clusters over the last 160 Myr with 
a total initial mass in the clusters of ~ 10^ Mq. This large mass is still only a small fraction 
of both the total mass of molecular gas (Gao et al. 2001) and the total luminous mass in the 
Antennae. Thus, the amount of gas available to fuel star formation in the Antennae does 
not require that the current epoch of intense star cluster formation be confined to a very 
short period of time. A model with a roughly constant cluster formation rate is consistent 
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with a picture in which the starburst activity in the Antennae is dynamically triggered, in 
which case we would not expect the star formation rate and star cluster formation rate to 
change on a much shorter timescale than the dynamical time, or on the order of 100 Myr 
(Mihos et al. 1993). 

5.2. The Future Rate of Cluster Formation in the Antennae 

We estimate the future rate of star cluster formation in the Antennae using a model 
based on our knowledge of current star formation in the Milky Way. If super giant molecular 
complexes are similar in their star formation properties to giant molecular clouds in the 
Milky Way, then we might assume that they will form at most 5% of their mass into stars. 
By analogy with what is seen in the Orion molecular cloud (Lada ct al. 1991a), wc would 
expect these stars to form predominantly in star clusters. Given the total mass in super 
giant molecular complexes is 6 x 10^ Mq, we might expect them to form a total mass of 
star clusters of up to 3 x 10^ M0. Zhang & Fall (1999) find that the mass function of the 
young massive clusters in the Antennae follows a power law with a slope of -2. This mass 
function implies that there is an equal amount of mass in equal logarithmic mass bins and 
thus the number of clusters more massive than 10^ Mq depends on the lower mass cutoff to 
the mass function. For example, if star clusters form over a mass range from 10^ to 10^ Mq, 
then the mass function will contain 600 clusters more massive than 10^ Mq, whereas if star 
clusters form over a mass range from 100 to 10^ Mq, the number of massive clusters would 
be only half as large. These rough estimates suggest that we might expect the existing super 
giant molecular complexes in the Antennae to form between 300 and 600 star clusters more 
massive than 10^ Mq. These clusters will likely form over at least a crossing time, which for 
the largest complex is 30 Myr. This imphes a massive star cluster formation rate of 10-20 
Myr"-*^, which is approximately the same rate as is seen in the youngest star cluster sample 
of Zhang & Fall (1999). 

If clusters older than 30 Myr have been subject to destruction, it is possible that the 
Antennae have been experiencing a high and uniform rate of star cluster formation over 
the last 100-200 Myr and could possibly sustain this formation rate for at least another 30 
Myr. What might happen after that? If the molecular gas can reform itself into super giant 
molecular complexes (or other structures) capable of forming massive star clusters over that 
same timescale, then massive star cluster formation could continue at the same rate until 
all the gas is exhausted. However, the north-east region of NGC 4038 suggests that, once a 
region has formed many young massive star clusters, it is unlikely to be ready to form them 
again just 30 Myr later. On the other hand, the fact that we see star clusters forming now 
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in the overlap region suggests a mode of star formation where a portion of the galaxy lights 
up in young clusters formed 30 Myr ago, while another portion of the galaxy is forming the 
next generation of clusters. Perhaps the western arc represents an intermediate age region 
where we still see gas and star clusters intermingled. In this picture, a continuous period of 
star cluster formation extending well in excess of 30 Myr but isolated to only a portion of 
the galaxy at a time could perhaps be sustained. If this mode of cluster formation continued 
for as long as 600 Myr, the Antennae would form 6000-12,000 clusters, of which as few as 
300-600 might survive in the long term. Slow infall of atomic material from the tidal tails 
(Hibbard et al. 2001; Hibbard & Mihos 1995) could help prolong the period of active star 
cluster formation. 

We need to consider whether the expected total star cluster formation rate in this 
simple model would cause the Antennae to exceed the typical specific frequency for globular 
(and other) clusters in an eUiptical galaxy. The current F-band absolute magnitude of the 
Antennae system is about -22 (from data in LEDA). The "normal" value for the specific 
frequency is about 3.5 (Harris 2001), which means we would expect to find a total of about 
2200 clusters more massive than 10^ Mq in an elliptical galaxy with My = —22 mag. If 
the Antennae were to fade by 1.5 mag in the process of forming an elliptical (Whitmore 
& Schweizer 1995), then the expected number of clusters would be 600. Whether these 
numbers agree in detail with the model discussed above depends upon both the duration 
of cluster formation and how many of the young massive star clusters survive over the long 
term. There is no difficulty if only 5-10% of the young massive star clusters that could be 
formed in future from the molecular gas reservoir survive as clusters after a Gyr. However, 
if all the young (2-6 Myr) massive star clusters that wc sec now will survive to old age, 
we would have to be witnessing an extremely short-duration starburst in order for the final 
specific frequency of the Antennae not to exceed (by a wide margin!) the typical specific 
frequency seen in elliptical galaxies. In this context, it is relevant that somewhat older (~ 1 
Gyr) mergers such as NGC 3921 (Schweizer et al. 1996) and NGC 7252 (Miller et al. 1997) 
are likely to end up with specific frequencies that fall in the normal range of field elliptical 
galaxies. 



6. Forming Massive Steir Clusters: Constraints from the Gas Properties 

6.1. Two Contrasting Models for Massive Star Cluster Formation 

Larsen & Richtler (2000) show that the fraction of the U-band luminosity that orig- 
inates in young star clusters is proportional to the star formation rate per unit area in a 
galaxy. Roughly speaking, their result suggests that if you double the surface density of 
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star formation, you double the fraction of the galaxy luminosity that occurs in star clusters. 
The luminosity in star clusters depends in a complicated way on the star formation history 
of a galaxy as well as on the survival rate of the star clusters. However, using the U-band 
luminosity function will preferentially select clusters that are typically < 30 Myr old, which 
is comparable to the timescale over which the star formation rate is traced via Ha imaging. 
Larsen (2000) concludes that galaxies form young massive star clusters wherever the star 
formation rate is high enough. He suggests that this could be hnked to the formation of 
Giant Molecular Associations, which may assemble more easily in galaxies with a higher gas 
surface density. However, our comparison of M51 and the Antennae in §4 shows that they 
have similar gas surface densities and yet have very different star cluster populations, and 
so gas surface density cannot be the only important parameter in the formation of young 
massive star clusters. 

There have been various models proposed to explain the formation of young massive 
star clusters. Schweizer et al. (1996) has noted the similarity between the masses of the 
young star clusters and the masses of giant molecular clouds in the Milky Way (10^ — 10^ 
Mq). However, the radii of the two objects are quite different, with the young massive star 
clusters having effective radii of 4 pc (Whitmore et al. 1999) and giant molecular clouds 
having radii in the range of 5 to 50 pc, depending on their mass. Schweizer et al. (1996) have 
suggested that a young massive star cluster could be formed through shock compression of 
a giant molecular cloud to force the cloud to form stars with an efficiency of ~ 50%, much 
higher than the star formation efficiencies of a few percent that are typically seen (Evans & 
Lada 1991). They suggest that one means by which giant molecular clouds could possibly 
be compressed is through the formation of a hot, high-pressure interstellar medium during 
galaxy collisions (Jog & Solomon 1992). This model was originally developed to explain 
the formation of a starburst during galaxy interactions. In the model, the outer layers of a 
giant molecular cloud arc compressed by contact with hot, shocked HI clouds and form stars 
with a high efficiency However, the amount of the cloud mass that is subject to this high 
efficiency is estimated to be 10%, and would be distributed in a shell-hke structure rather 
than a compact, spherical structure. In addition. Jog & Solomon (1992) assume that only 
stars more massive than 1 Mq are formed in this process in order to match the observed 
parameters of starburst galaxies, in particular their high infrared luminosity. Thus, this 
model does not seem to be directly applicable to the formation of compact massive star 
clusters. 

An alternative model that has been proposed to explain the formation of globular clus- 
ters (Harris & Pudritz 1994; McLaughlin & Pudritz 1996) may be applicable to the formation 
of young massive star clusters in galaxy mergers. This model supposes that globular clusters 
formed by the same type of process that we see forming star clusters today. However, the 



-20- 



large masses of globular clusters compared to star clusters in clouds like Orion requires that 
the globular clusters form within molecular clouds that are substantially larger than the 
molecular clouds in the Milky Way. For example, the Orion B molecular cloud has a total 
mass of 8 X 10*^ Mq and contains five molecular cores with masses ranging from 100 to 450 
Mq (Lada, Bally, & Stark 1991b). Four of these cores are forming infrared star clusters 
with a star formation efficiency of 50% (Lada et al. 1991a) and most of the star formation is 
occurring in these cores. If a globular cluster (or a young massive star cluster) formed inside 
a larger molecular core with a star formation efficiency of 50%, this would require cores with 
masses from 2 x 10^ to 2 x 10^ Mq. If these large cores were contained within a larger cloud 
with a core-to-cloud mass ratio similar to Orion B, the cloud mass required to contain these 
cores would be 4 x 10^ to 4 x 10^ M©. These cloud masses are very similar to those of the 
super giant molecular complexes we have identified in the Antennae. Note that in this model 
we would expect each super giant molecular complex to contain several massive cores and 
hence to form several massive star clusters over its lifetime. 



6.2. The Role of Pressure in the Two Models 

One possible way to distinguish between these two models may be via the environment 
in which the clouds live, particularly the pressure. We can estimate the pressure in a few 
regions in the Antennae from the recent work by Fabbiano et al. (2003). They have made 
Chandra observations which are sufficiently sensitive to be able to fit the physical conditions 
in the hot interstellar medium in the two nuclei and also in a region in the western star- 
forming arc. Assuming a scale height of 200 pc for the high temperature component, the 
pressure in each of these regions is in the range of 4 — 8 x 10^ K cm~^, with an average over 
the three regions of 6 x 10^ K cm~^. This pressure is significantly larger than the typical 
interstellar medium pressure in the Milky Way of 10^ K cm~^ (Elmegreen 1989). If the 
pressure in these three regions is typical of the Antennae disk as a whole, then the super 
giant molecular complexes in the Antennae are embedded in a higher pressure environment 
than are molecular clouds in the Milky Way. 

How does this observed pressure compare to the models of star cluster formation dis- 
cussed previously? Harris & Pudritz (1994) adopt a surface pressure for molecular clouds 
that is ten times larger than the typical interstellar medium pressure to account for the 
effect of embedding HI envelopes. Using the formula in Elmegreen (1989) and HI data from 
Hibbard et al. (2001), the pressure due to the average atomic gas column density in the An- 
tennae is only ~ 10^ K cm~^, closer to the Milky Way pressure than any other component of 
the interstellar medium in the Antennae. Thus, the surface pressure felt by the super giant 
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molecular complexes in the Antennae is likely to be in the range of 6 x 10^ to 6 x 10^ K 
cm~^, depending on whether they are surrounded by massive HI envelopes with pressures 
significantly above the average pressure in the atomic gas. Since the radius of a cloud varies 
as P"^/'^, a cloud of 3 x 10^ Mq should have a radius between 300 and 600 pc. This size is a 
reasonable match to those of the molecular complexes in the Antennae. The model of Jog & 
Solomon (1992) predicts a pressure of 10® K cm~^ in the hot, shocked HI clouds. Although 
this pressure is substantially higher than that estimated from the X-ray data, the expected 
filling factor of this hot gas is less than 1%. With this low filling factor, the average pressure 
predicted in the model is 6 x 10^ K cm~^, which is quite similar to the X-ray data. Thus, 
it seems that the pressure of the interstellar medium in the Antennae cannot help us to 
distinguish between the two models. 

Even the high pressures inferred from the X-ray data for the Antennae are still substan- 
tially smaller than the typical pressure inside a single, self-gravitating molecular complex. 
The average pressure in SGMC-3 is 2 x 10^ K cm~^, while the central pressure is in the 
range of 10^ — 10^ K cm~^, depending on the exact geometry. Indeed, the pressure inside 
the p Ophiuchus core, a relatively small star-forming cloud in the Milky Way, is 10^ — 10® K 
cm~^ (Johnstone et al. 2000). Elmegreen & Efremov (1997) estimate the pressure required 
to form a typical globular cluster to be in the range of 10® — 10® K cm~^. 



6.3. Applying the Models to the Non-Merger Environment 

It is important to see if we can understand the formation of young massive star clusters 
in non-interacting galaxies in the context of these two models. A particularly useful case 
study is the R136 star cluster in the LMC The total mass of stars in R136 is estimated to be 
6 X 10"^ Mq (Hunter et al. 1995) and the age of the star cluster is at most 1-2 Myr (Massey 
& Hunter 1998). There appears to be no molecular gas in the immediate environment of 
R136 (Johansson et al. 1998); however, there is a very extended region of CO emission to 
the south of R136 (Cohen et al. 1988). If we consider this entire region as a single large 
structure, it has a size of 400 x 1200 pc and a total mass of about 9 x 10^ Mq (adopting the 
appropriate CO-to-H2 conversion factor from Wilson (1995)). Adopting the velocity width 
of 10.5 km s~^ measured by Kutncr et al. (1997), the virial mass of this region is 1.5 x 10^ 
Mq or within a factor of two of the flux-based mass. Thus, this large CO complex may be 
gravitationally bound. If R136 formed in a core of 1 x 10^ Mq, wc would expect this core 
to have lived inside a larger cloud of mass 2 x 10^ Mq, which is similar to the mass of the 
large CO complex. Thus, it seems reasonable to suggest that R136 might have formed from 
the most massive core of that molecular complex. R136 and 30 Doradus are rather isolated 
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at present from the main part of this molecular complex. However, the intervening region 
contains a number of HII regions and supernova remnants (Cohen et al. 1988), which could 
have acted along with R136 and 30 Doradus to destroy much of the molecular gas over this 
larger region. Thus, the environment of the R136 cluster appears roughly consistent with the 
cluster formation picture put forward by Harris & Pudritz (1994), while the over-pressure 
model advanced by Schweizer et al. (1996) seems less applicable here. 

Returning now to M51, the CO data show that it currently contains 13 molecular 
complexes with masses in the range of 2 — 5 x 10^ Mq (Rand & Kulkarni 1990). Scaling 
again from Orion, we would expect the most massive core in each complex to be in the 
range of 10^ — 3 x 10^ Mq and hence to be able to form a young star cluster with a mass 
up to 1 - 2 X 10^ Mq. The 69 star clusters identified by Larsen (2000) have ages <500 
Myr. Depending on the masses of these star clusters (which have not been estimated) it 
seems quite possible for this relatively small population of clusters to have been formed from 
objects like the lower mass molecular complexes seen in the Antennae. 

Young massive star clusters have been found even in relatively low-luminosity dwarf 
irregular galaxies (Gelatt, Hunter, & Gallagher 2001; Hunter et al. 2000; Billett, Hunter & 
Elmegreen 2002). Billett et al. (2002) note that, although it is rare for a dwarf galaxy to 
form luminous star clusters, when they do, they tend to form several clusters with similar 
ages. This result suggests that star cluster formation is concentrated to a localized region. 
Billett et al. (2002) emphasize that a lack of shear is probably the most important difference 
between dwarf irregular and spiral galaxies, and suggest that triggered large-scale flows, 
possibly by an interaction, or ambient instabilities in a shear-free environment can make the 
clouds that form young massive star clusters in dwarf galaxies. Given the radically disturbed 
velocity fields seen in the Antennae, reduced shear might also be a factor in cloud formation 
in this system. An important difference between the Antennae and dwarf galaxies is the 
sheer size of the gas reservoir (Gao et al. 2001); perhaps combining a large gas reservoir 
with localized regions of reduced shear is the key to explaining the explosion of star cluster 
formation within the Antennae. 



7. Conclusions 

We have used sensitive CO J=l-0 observations of the Antennae to study the detailed 
properties of the molecular clouds with the goal of understanding the prodigious formation 
of young massive star clusters in this nearby merger system. We have identified a total of 
~ 100 clouds in the data cube with masses ranging from 2 x 10® to 9 x 10^ Mq. This sample 
of extragalactic molecular clouds is unique in the total number of clouds identified and in 
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the mass range probed by the observations. 

Above our 5a completeness limit of 5 x 10^ Mq, the cloud mass function has a slope of 
— 1.4 ± 0.1. This mass function slope is very similar to that seen in molecular clouds and 
molecular cloud cores in the Galaxy, and is somewhat steeper than the mass function slope 
of -2 estimated for the luminous young star clusters in the Antennae by Zhang & Fall (1999). 
Our data suggest that the molecular interstellar medium in galaxies is governed by the same 
mass function slope over 8 to 9 orders of magnitude in mass, from 1-10 Mq up to 10^ Mq. 

We have compared the Antennae with the nearby spiral galaxy M51, for which similar 
sensitivity and spatial resolution CO observations exist. Although the two galaxies have 
similar gas surface densities and total gas masses, the molecular clouds in M51 are an order 
of magnitude less massive than those in the Antennae. In addition, M51 has a much smaller 
population of young massive star clusters, perhaps 50 times smaller than that of the An- 
tennae. One significant difference between the gas in the two galaxies is that M51 exhibits 
a smooth velocity field, while the velocity field in the Antennae is highly disordered. One 
possibility is that young massive star cluster formation in M51 is suppressed by its inability 
to form very massive gas clouds due to high shear in its differentially rotating disk. 

Comparing the CO data with the locations of the youngest clusters shows that many 
of these clusters lie in the CO-rich overlap region. The extremely high extinction {Ay ~ 100 
mag) towards the CO peaks means that the star clusters must he in front of most of the 
molecular gas in this region. Interestingly, some of these youngest clusters are found as much 
as 2 kpc from regions with detectable molecular gas. This result implies that either some 
young massive star clusters can form from clouds less massive than 5 x 10^ M©, or these 
clusters have already destroyed their parent molecular clouds, or both. 

We have combined our CO data with pubhshed radio and optical data to sketch out the 
recent star formation history of the Antennae and to speculate on how star formation will 
proceed in the future. The relative numbers of very young massive star clusters and thermal 
and non-thermal radio continuum sources are consistent with a constant star formation rate 
over the last 10-20 Myr. If the star formation rate has been constant over the last 100 Myr, 
then the relative numbers of star clusters with ages less than 10 Myr and 20-100 Myr implies 
that a very large fraction of the star clusters formed must evaporate or be destroyed (Zhang 
& Fall 1999). On the other hand, if most star clusters survive for long periods, we would 
have to be witnessing an extremely short-duration starburst in order for the final specific 
frequency of the Antennae not to exceed the typical specific frequency seen in elliptical 
galaxies by a wide margin. 

The abundant supply of molecular gas seen in the extremely massive molecular clouds 
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suggests that star cluster formation could easily proceed for at least a crossing time of one of 
these large clouds, or about 30 Myr. Indeed, the large reservoir of molecular gas measured 
by Gao et al. (2001) suggests that star cluster formation could continue for quite some 
time. The current distribution of gas and stars in the Antennae suggests a mode of star 
cluster formation where portions of the galactic disks are lit up with recently formed clusters 
(for example, regions like the north-western arc) while other regions are forming the next 
generation of star clusters (currently in the overlap region). 

We have compared our observations with two different models to explain the formation 
of young massive star chistcrs. The identification of extremely massive gas clouds in the 
Antennae means that the model of Harris & Pudritz (1994), which envisions globular cluster 
or super star cluster formation as a scaled up version of Galactic star cluster formation, is a 
potentially viable model. The model advanced by Schweizer et al. (1996), which requires a 
source of high-pressure in order to collapse pre-existing giant molecular clouds to form stars 
with a much higher efficiency, may also be viable given recent constraints on the pressure in 
the Antennae from X-ray observations (Fabbiano et al. 2003). In addition, higher pressures 
can also be a feature of the Harris & Pudritz (1994) model, so that pressure alone cannot 
distinguish between the models. Understanding the formation of super star clusters in dwarf 
galaxies, where there may not be a global source of enhanced pressure, may provide additional 
interesting constraints on this question. Ultimately, observations with new telescopes such 
as the Atacama Large Millimeter Array, which will have sufficient resolution and sensitivity 
to probe the small physical scales on which young massive star clusters form, will probably 
be required to fully understand this intriguing mode of star formation. 
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Institute of Technology and is supported by NSF grant AST96-13717. CDW acknowledges 
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A. Properties of Individual Super Giant Molecular Complexes in the 

Antennae 

The properties of the individual super giant molecular complexes (such as position, flux, 
and velocity) are given in Tabic 2. The clouds were identified using the program CLFIND 
with a contour level of 0.11 Jy beam^^ {2a) and their properties were measured using the 
program CLSTATS (Williams et al 1994). The 5a mass sensitivity limit of the data set is 
5 X 10^ M„. 



B. The CO-to-H2 Conversion Factor in the Antennae 

Most of the clouds identified in the Antennae are unresolved, and so their masses can 
only be calculated from their CO fiux. It is important, therefore, to estimate the CO-to- 
H2 conversion factor that is appropriate for the Antennae using the few clouds that arc 
large enough that their true diameters can be deconvolved from the synthesized beam. Wc 
excluded clouds in the two galactic nuclei or near bright mid-infrared peaks from this analysis 
because the conversion factor could be modified in these regions due to intense star formation 
or high pressure (Bryant & Scoville 1996; Solomon et al. 1997). 

Wc analyzed the data cube with CLFIND and CLSTATS using three difi^erent contour 
levels (0.10, 0.11=2cr, and 0.12 Jy beam^^) as described above. The integrated intensity map 
for each cloud that appeared resolved in both dimensions from the output from CLSTATS 
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was inspected to see if it was in fact resolved. (Clouds may appear resolved from their sizes 
but not actually be resolved if, for example, the cloud is actually made up of two peaks, or 
if the cloud is highly elongated in one direction.) We identified four, three, and two clouds 
that appeared to be resolved when identified with contour levels of 0.10, 0.11, and 0.12 Jy 
beam^^, respectively (Table 3). Only one cloud was found in common to both the 0.10 and 
0.11 Jy beam~^ processing. The remaining seven clouds were identified when contoured at 
different levels, but, due to small changes in how flux was assigned to each cloud, they did 
not always appear to be resolved spatially. 

For each resolved cloud, the radius and velocity width were measured using CLSTATS 
while the total CO flux was measured from a gain-corrected zeroth moment map of each 
individual cloud. We used these measurements to calculate both the virial mass and the 
molecular mass. The virial mass is given by 

where ^Vfwhm is the full width half maximum velocity in km s~^ and Rpc — ^ Areaj'K is 
the deconvolved radius of the cloud in parsecs. The molecular mass is given by 

M^o, = 1.61 X 10^Dl,p,5co Mo 

where Dupc is the distance to the cloud in Mpc and Sco is the CO integrated intensity in Jy 
km s~^ (Wilson & Scoville 1990). The ratio of the assumed value for CO-to-H2 conversion 
factor in the Galaxy relative to the true value in the Antennae can then be determined from 
the ratio of the molecular mass to the virial mass. The average value of this ratio for the 
eight resolved clouds from the three different identification runs is 1.3 ±0.3 with a standard 
deviation of 0.7. This result suggests that the CO-to-H2 conversion factor may be slightly 
smaller in the Antennae than in the Milky Way. However, given that the ratio is consistent 
with equal conversion factors within the uncertainties, we have chosen to use the standard 
Galactic value (3 x 10^° H2 cm ^ (K km s ^) ^) for ease of comparison with other work. 

Clouds that were resolved but were not used to determine the CO-to-H2 conversion 
factor due to their location in the nuclei or the mid-infrared bright part of the overlap region 
are given in Table 4. Cross-identifications with Paper I are also given; note that, in some 
cases, several complexes in Table 2 and Table 4 correspond to a single super giant molecular 
complex from Paper I. A comparison of the virial and molecular masses in Table 4 suggests 
that the CO-to-H2 factor may be smaller than the adopted value in several of these clouds. In 
particular, the complexes in the nucleus of NGC 4038 all have molecular masses that exceed 
their virial masses by at least a factor of three. This result suggests that the CO emission 
is over-luminous in this galactic nucleus, perhaps similar to the effect seen in ultraluminous 



-31 - 



infrared galaxies (Solomon et al. 1997). We will discuss the properties of the nuclear regions 
in a future paper (Wilson, Madden, & Charmandaris, in preparation). 
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Table 1. Massive Molecular Complexes in Other Galaxies 



Galaxy 


Beam Size 


Distance 


Maximum Mass 


Notes 


References" 




(pc) 


(Mpc) 


(108 Mo) 






M83 


110 X 55 


3.3 


0.1 


10 complexes 


1,10 


NGC 5055 


220 X 170 


7.2 


0.5 


17 complexes 


2,11 


M51 


370 X 290 


8.6 


0.5 


26 complexes 


3,12 


NGC 1068 


210 X 210 


14.4 


4.5 


38 complexes 


4,11 


MlOO 


360 X 280 


16 


1.2 


7 complexes 


5,13 


NGC 4414 


480 X 480 


19 


0.4 


8 complexes 


6,14 


NGC 4038/39 


480 X 310 


19 


8.9 


86 complexes 


7,15 


Arp 299 


750 X 510 


42 


40 


nucleus of IC694 


8,15 


Arp 220 


210 X 190 


77 


10 


nucleus, Xco uncertain 


9,15 



Note. — Masses have been scaled to the same CO-to-H2 conversion factor used in this paper, except for 
Arp 220 where a value of 2.5 times smaller was used. 

"References for molecular data and distance. 1. Rand et al. (1999) 2. Thornley & Mundy (1997) 3. 
Rand & Kulkarni (1990) 4. Planesas et al. (1991) 5. Rand (1995) 6. Sakamoto (1996) 7. this paper 8. 
Casoli et al. (1999) 9. Sakamoto et al. (1999) 10. Distance to group member NGC 5253 from Gibson ct 
al. (2000) 11. TuUy & Fisher (1988) 12. Feldmeier et al. (1997) 13. Ferrarese et al. (1996) 14. Turner et 
al. (1998) 15. Ho = 75 km s"! Mpc"!. 
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Table 2. Giant Molecular Complexes in the Antennae 



ID*' 


a 


5 


Sco 




Vlsr 


AVpwHM 


Mmol 




(2000) 


(2000) 


(Jy km s~^) 


(K) 


(km s ^) 


(km s ^) 


(10^ Mq) 


20 




-18'53'1 5 

XO.tJfJ. XiJ 


16.2 


3.1 


1683 


21 


9.4 


86 


1 9-01 ■'^'^ 1 


- XO.OO. XU 


1.8 


2.6 


1755 


Q 


1.0 


22° 


X Zr . X • ij (J • J. 


-18:53:14 


23.2 


4.0 


1703 


37 


13.5 


58° 


1 9-01 9 


xo.uo. xo 


47.4 


2.8 


1610 


56 


27.5 


72 




XO.OO. XU 


0.4 


2.1 


1641 


5 


0.2 


32 




XO.tJO.lJO 


2.9 


1.9 


1464 


16 


1.7 


79 


12:01:53.7 


-18:53:08 


5.5 


2.0 


1683 


27 


3.2 


65 


1 2-01 -53 7 


-18:53:07 


0.6 


1.6 


1620 


9 


0.3 


7° 


1 2-01 -54 9 


-1 S-53-n6 


11.2 


3.5 


1501 


11 


6.5 


80 




-1 K-'^'^-D'^ 
-* xo.oo.uu 


0.9 


1.9 


1693 


9 


0.5 


44 


1 2'ni '52 4 


_i «.c;Q.n4 


1.4 


2.5 


1558 


7 


0.8 


69 




-18-53'n4 

XO.iJO.U^ 


1.3 


2.8 


1631 




0.7 


13 




xo.oo.uo 


2.1 


2.3 


1480 


5 


1.2 


12 


1 9-ni ■'^4 fi 


-18:53:02 


40.6 


2.8 


1485 


53 


23.6 




1 2-01 -55 


-18:53:01 


42.6 


3.8 


1532 


31 


24.8 


15° 


12:01:54.9 


-18:53:01 


22.7 


2.6 


1568 


49 


13.2 


9° 


1 2-01 -54 9 


-IS-^'^-ni 

iO.tJiJ.VJX 


11.3 


3.7 


1553 


12 


6.6 


66 




xo.uz.oy 


4.7 


1.8 


1620 


22 


2.7 


34 


1 2'ni '55 4 


-18* ^2* ^8 

XO<i-f Z.<_IO 


5.4 


2.2 


1537 


15 


3.1 


38 


1 2'01 '55 8 

X Zi . X . ij . o 


-18'52'58 


0.4 


2.3 


1542 


g 


0.2 


77 


XZi. L/X . Ut:. I 


- XO .tJZ . 0\J 


2.5 


1.5 


1667 


12 


1.5 


54 


1 9-01 -^d Q 

X Zi . L/ X . ij t: . C/ 


- XO .OZ .tJU 


1.4 


2.0 


1584 




0.8 


63 


1 9-ni -^d Q 


1 S■'^9■f^'^ 
- xo.oz .ou 


0.4 


1.5 


1605 




0.2 


73 


1 2'01 '54 7 


XO.tJZ.tJiJ 


1.8 


1.5 


1641 


18 


1.0 


14 


1 2'01 '54 8 


-18'52'54 


4.3 


2.1 


1521 


55 


2.5 


40 


X Z ■ U X . OtJ. X 


XO . tJZ . tj'rt 


0.3 


1.6 


1547 


5 


0.2 


33 


1 2-01 -55 4 

XZ»UX»Ot^»rr 


-1 S-52-54 


1.6 


1.9 


1475 


5 


0.9 


59 


1 2-01 -54 8 

XZ»L/X»v.Jrr»0 


-1 S-52-53 

X (J . tj Z . I J 


2.4 


2.0 


1594 


17 


1.4 


70 


1 9'ni •'=i4 7 

XZ.Lf X . 1 


-1 K-f^9-f^'^ 


0.4 


1.4 


1631 


Q 


0.2 


4° 


1 2'ni '54 8 

X Zi . U X ■ tJT: ■ O 


-1S*^9*^9 

XO*i_fZ*i_IZ 


23.8 


2.8 


1449 


23 


13.8 


5° 


1 9-01 '^d Q 


XO.tJZ.tJZ 


44.0 


3.4 


1480 


33 


25.6 


47 


1 9*01 '"^d 7 

X Z a X • ^Jrr • 1 


-18'52'52 

XO»tJZ»OZ 


2.3 


1.4 


1563 


13 


1.3 


O / 


Iz.Ul.O 1. ( 


- lo. Dz. Dz 


n o 
u.y 


1 c: 
i.o 


iuo / 




n c: 


11 


12:01:54.7 


-18:52:50 


8.6 


2.2 


1417 


23 


5.0 


6 


12:01:55.4 


-18:52:49 


109.0 


4.1 


1501 


39 


63.3 


25 


12:01:54.8 


-18:52:49 


3.5 


1.6 


1381 


13 


2.1 


26 


12:01:55.4 


-18:52:49 


1.4 


1.8 


1391 


16 


0.8 


27 


12:01:54.7 


-18:52:49 


1.4 


1.4 


1397 


15 


0.8 


29 


12:01:55.4 


-18:52:48 


6.1 


1.9 


1428 


15 


3.6 


31 


12:01:55.4 


-18:52:48 


3.6 


2.0 


1459 


8 


2.1 


24 


12:01:55.4 


-18:52:48 


1.5 


2.0 


1381 


8 


0.9 


23 


12:01:55.4 


-18:52:47 


8.3 


2.3 


1371 


22 


4.8 


30 


12:01:55.4 


-18:52:47 


2.2 


2.0 


1443 


11 


1.3 


83 


12:01:55.4 


-18:52:47 


1.4 


1.8 


1709 


6 


0.8 


28 


12:01:55.5 


-18:52:46 


3.4 


2.0 


1402 


10 


2.0 


10 


12:01:54.7 


-18:52:44 


14.9 


2.8 


1589 


23 


8.6 


67 


12:01:54.4 


-18:52:44 


5.7 


1.8 


1620 


13 


3.3 


74 


12:01:54.5 


-18:52:44 


0.6 


1.4 


1646 


5 


0.3 
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Table 2 — Continued 



ID'' 


a 


S 


Sco 






^Vfwhm 






(2000) 


(2000) 


(Jy km s~^) 


(K) 


/I — 1 \ 
(km s ^) 


/I — 1 \ 
(km s ^) 


(10^ M0) 


16 




-18:52:41 


16.0 


2.5 


1568 


21 


9.3 


17 




-1 8-52-40 


11.4 


3.1 


1584 


14 


6.6 


18 


1 2'ni •'i5 1 


-1 8-52-38 


7.3 


2.9 


1599 


13 


4.3 


19 


1 9-ni ■'^4 7 


XO. UZi. 


23 6 


2.2 


1620 


26 


13.7 


35 


1 9-ni -^"^ 9 

J. ■ U -L . o o . ^ 


1 K-'^9-'^'^ 

- LO'Ozi, OO 


0.6 


1.5 


1527 


5 


0.4 


84 


1 2-ni -^5 


-18:52:33 


1.9 


1.8 


1719 


g 


1.1 


53 


1 2-01 -53 




1.6 


1.7 


1579 


7 


0.9 


76 




-1 8-52-31 


2.1 


2.5 


1651 


Q 


1.2 


75 


1 2'ni •'n4 4 


-1 8-52-30 

XO. tJZl.tJKJ 


2.0 


1.4 


1646 


17 


1.2 


51 


12:01:54.4 


-18:52:29 


9.2 


1.6 


1568 


44 


5.4 


55 


1 2-ni -55 1 


-18-52-29 


4.7 


2.1 


1589 


7 


2.7 


42 


1 2-ni -53 5 


-18-52-28 


0.8 


1.4 


1553 


5 


0.5 


68 


12:01:54.5 


-18:52:28 


3.5 


1.5 


1625 


16 


2.1 


52 


1 2'ni '55 1 


-18-52-26 


2.2 


1.9 


1568 




1.3 


60 


12-01 -55 9 


-18:52:25 


1.7 


2.4 


1594 


g 


1.0 


71 


1 2'01 '50 5 


-1 8-52-24 


16.5 


1.9 


1667 


41 


9.6 


45 


12-01 -51 Q 


-1 8-52-22 


0.7 


1.6 


1558 


5 


0.4 


78 


12:01:52.0 


-18:52:20 


4.1 


1.6 


1677 


13 


2.4 


56 


12-01 -53 3 


-18-52-16 

XO»iJ^» xu 


2.1 


1.4 


1584 


22 


1.2 


48 


1 2-01 -"13 7 


-1 8-52-1 2 

~XO-^^ - x^ 


1.3 


1.4 


1563 


g 


0.7 


64 


1 2-01 -"12 1 


-1 8-52-08 

XO. tJZi.\J(j 


4.7 


1.3 


1615 


28 


2.7 


62 


1 2-01-50 Q 


-1 8-52-07 

XO. 1 


5.4 


1.6 


1599 


19 


3.2 


57 


12-01-50 8 


-18-52-06 


2.4 


1.4 


1589 


12 


1.4 


61 


12:01:50.5 


-18:52:06 


0.6 


1.5 


1594 


5 


0.4 


49 


12:01:53.1 


-18:52:05 


1.5 


1.6 


1563 


8 


0.9 


81 


12:01:53.8 


-18:52:04 


1.7 


1.4 


1693 


9 


1.0 


1 


12:01:53.0 


-18:52:02 


152.4 


6.7 


1620 


45 


88.6 


2 


12:01:53.0 


-18:52:01 


50.8 


5.8 


1651 


16 


29.5 


3 


12:01:53.0 


-18:52:01 


35.9 


4.5 


1667 


18 


20.9 


21 


12:01:53.0 


-18:52:00 


19.5 


2.4 


1688 


27 


11.3 


50 


12:01:51.8 


-18:51:58 


1.2 


1.3 


1563 


8 


0.7 


39 


12:01:50.7 


-18:51:57 


9.1 


1.7 


1547 


39 


5.3 


82 


12:01:52.3 


-18:51:57 


4.7 


1.5 


1703 


20 


2.7 


36 


12:01:50.7 


-18:51:54 


7.1 


1.5 


1532 


25 


4.1 


85 


12:01:52.3 


-18:51:54 


1.9 


1.6 


1729 


7 


1.1 


41 


12:01:51.1 


-18:51:52 


17.2 


2.1 


1568 


29 


10.0 


46 


12:01:51.8 


-18:51:41 


4.5 


2.3 


1558 


11 


2.6 


43 


12:01:51.8 


-18:51:38 


4.9 


2.5 


1553 


17 


2.8 



Note. — A distance to the Antennae of 19 Mpc is assumed throughout. 
^Cloud ID number assigned by CLFIND algorithm. 

''Observed peak brightness temperature excess above the 2.74 K cosmic background. 

'^Tlicsc clouds belong to the super giant molecular complexes discussed in Paper I: complexes 4, 
5: SGMC 2; complexes 7, 8, 9: SGMC 4; complex 15: SGMC 5; complexes 22, 58: part of NGC 
4039. 
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Table 3. Resolved Clouds Used to Determine Xqo 



ID 


Sco 


R 


AVpwHM 






CLFIND contour 




(Jy km s~^) 


(pc) 


(km s"-"-) 


(108 Mq) 


(108 Mq) 


(Jy beam^^) 


17 


28.6 


580 


33 


1.3 


1.7 


0.12 


39 


9.6 


430 


23 


0.45 


0.56 


0.12 




23.6 


560 


24 


0.64 


1.4 


0.11 


36 


7.1 


430 


25 


0.53 


0.41 


0.11 


71 


16.5 


610 


41 


2.0 


0.96 


0.11 


41 


16.9 


610 


20 


0.48 


0.98 


0.10 


62 


19.3 


600 


46 


2.5 


1.1 


0.10 


46 


5.7 


410 


14 


0.16 


0.33 


0.10 



Note. — A distance to the Antennae of 19 Mpc is assumed throughout. Values for Sco and 
AVfw H M may differ from Table 2 for the clouds identified with a different contour level. 

^This cloud was resolved for both the 0.11 and 0.10 Jy beam~^ contour levels. The values 
for all parameters are the average of two values. 



Table 4. Other Resolved Clouds 



ID 


R 


AVfwhm 






ID from Paper I 




(pc) 


(km s-i) 


(108 Mq) 


(108 Mo) 




1 


760 


45 


3.1 


8.9 


NGC 4038 


2 


600 


16 


0.3 


3.0 


NGC 4038 


3 


630 


18 


0.4 


2.1 


NGC 4038 


6 


670 


39 


2.0 


6.3 


SGMC 1 


12 


670 


53 


3.7 


2.4 


SGMC 3 


20 


440 


21 


0.4 


0.9 


NGC 4039 


21 


480 


27 


0.7 


1.1 


NGC 4038 
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Fig. 1. — (a) Integrated CO intensity map from Paper I. Contour levels are 
1,2,4,6,8,10,15,20,25.. .50,60.. .,90 Jy beam"! km s-\ This map was made using CLPLOT 
and including only molecular complexes with velocity widths greater than or equal to 15.6 
km s~^ found with CLFIND. This image has been corrected for the fall-off in sensitivity 
due to the primary beam, (b) Integrated intensity map made using the same data set and 
criteria, but including the clouds that were missed due to a bug in the CLFIND program, (c) 
A standard zeroth moment map made from the same data cube used to identify molecular 
complexes. The 50% gain limit of the telescope primary beam is clearly visible in this figure. 
Negative bowls exist in this map (see text) but are not shown for ease of comparison with 
the other two maps. 

Fig. 2. — First moment map of the CO emission in the Antennae. The color scale runs from 
1460 km s~^ to 1660 km s~^, while the contour levels run from 1460 to 1660 km s~^ in steps 
of 20 km s~^. Note in particular the disturbed velocity field in the overlap region. 

Fig. 3. — The differential mass function for the molecular complexes in the Antennae. Mass 
functions arc shown for three identification runs with different contour levels; each identifi- 
cation run is also shown for two different binnings of the mass function separated by 0.1 in 
log(M). The 5(7 completeness limit is indicated by the dashed fine. The contour level in Jy 
beam~^ and the slope derived from a least squares fit to the points above the completeness 
hmit are given in each panel. The average slope of —1.4 ± 0.1 is very similar to the slope 
observed for Giant Molecular Clouds (10^ - 10^ M©) in the Milky Way. 

Fig. 4. — Locations of young, red super star clusters using the coordinates from Zhang et al. 
(2001) overlaid on the CO contours from Figure lb. Roughly 40% of these clusters are found 
in the overlap region, which covers only 10% of the area of the total CO map. However, some 
young clusters are found as far as 1-2 kpc from strong CO emission, which suggests they 
either formed from a lower mass molecular complex or they have already destroyed their 
parent cloud. The locations of the eight resolved clouds from Table 3 are also indicated. 



